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(54) METHOD FOR DETECTING A PLURALITY OF SYMBOLS 

(57)Abstract: 

PROBLEM TO BE SOLVED: To provide the simple type of a 
detection method by a sphere in a maximum likelihood multi-user 
detecting method for a DS-CDMA(Direct Sequence Code Division 
Multiple Access) system. 

SOLUTION: In the method for detecting a plurality of symbols, each 
symbol belongs to a modulation constellation and is the objects of 
spectral spreading by means of a spreading sequence. The method 
comprises a filtering step for supplying a real vector characteristic 
of the received signal; the closest neighbor of the vector is sought 
within the lattice of points generated by the modulation 
constellation; the search is limited to candidate vectors and each of 
whose components has a value lying within a search interval defined 
by a lower bound and an upper bound; and the bounds are selected, 
so that each of the intervals comprises only component values of 
points situated within a sphere with a predetermined radius and 
belonging to the modulation constellation. 
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1. Title of the Invention 

METHOD OF DETECTING A PLURALITY OF SYMBOLS 

2. Claims 

1. Method of detecting a plurality of symbols (d k (/)) transmitted by or for a 
plurality K of users, each symbol belonging to a modulation constellation and being 
the subject of a spectral spreading by means of a spreading sequence, the said method 
comprising a filtering step adapted for supplying a real vector (z) characteristic of the 
said received signal, characterised in that at least the closest neighbour of the said 
vector is sought within a lattice of points (S) generated by the said modulation 
constellations, the said search being limited to candidate vectors (x), each of whose 
components (b t ) has a value lying within a search interval defined by a lower bound 
(BT ) and an upper bound the said bounds being chosen so that each of the said 
intervals comprises only component values of points situated within a sphere with a 
predetermined radius (JC) and belonging to a modulation constellation. 

2. Detection method according to Claim I, characterised in that, for a 
component with a given index f, each pair of bounds, lower (BT ) and upper (B?\ is 
calculated from a quantity T t characteristic of the width of the search interval for the 
said component, this characteristic quantity being determined by recurrence on the 
index fc the characteristic quantity (73 for a given index being determined from the 
quantity with the previous index and the value of a component chosen in 
the search interval ([BT.x&i] ) relating to the component with the previous index. 

3. Detection method according to Claim 2, characterised in that, for a candidate 
vector component with a given index (2k\ the lower bound of the search interval is 
chosen so as to be greater than the smallest integer {Mu ) corresponding to a symbol of 
the modulation constellation (A A ) and the upper bound of the search interval is chosen 
so as to be less than the largest integer (Afa) corresponding to a symbol of the 
modulation constellation. 

4. Detection method according to Claim 2, characterised in that, for a first 
candidate vector component with a given index the lower bound of the search 
interval is chosen so as to be greater than a first integer (Mit-i) and the upper bound of 
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the search interval is chosen so as to be less than a second integer (Jl/i-i), the said first 
and second integers being determined from the value of a second component (M of 
the candidate vector such that the first and second components relate to one and the 
same user. 

5. Detection method according to Claim 4, characterised in that the said first and 
second integers are respectively determined as the smallest value and the largest value 
of the first component such that the complex number defined by the value of the said 
first component and the value of the said second component is a symbol of the 
modulation constellation of the said user. 

6. Detection method according to one of the preceding claims, characterised in 
that the search for the closest neighbour is effected by scanning one by one the said 
search intervals for the different component indices (/) and choosing a component 
value (b,) in each of the said intervals, the bounds of each interval being determined 
according to the width of the interval with the preceding index and the 
component value (b m ) chosen in this same interval. 

7. Detection method according to one of the preceding claims, characterised in 
that, if, during the said search, the norm (|wj) of a candidate vector is less than the 
radius of the said sphere, the radius is updated to the value of the said norm 

3. Detailed Description of Invention 

The present invention concerns a multiuser detection method and device. More 
particularly, the present invention concerns a maximum likelihood multiuser detection 
method and device for a DS-CDMA (Direct Sequence Code Division Multiple Access) 
telecommunication system. 

In a DS-CDMA mobile telecommunication system, the separation of the 
communications coming from or going to the different users is achieved by multiplying 
each complex symbol of a user by a spreading sequence which is peculiar to him, also 
referred to for this reason as the user signature. The spreading frequency (chip rate) 
being greater than the frequency of the symbols, the signal transmitted by each user is 
distributed (or spread) in the space of the frequencies. The ratio between the band 
occupied by the spread signal and the band occupied by the information signal is 
referred to as the spreading factor. On reception, the separation of a given user is 
obtained by means of a filtering adapted to the corresponding signature. When the 
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transmission channel has a plurality of propagation paths, the output of the adapted 
filtering contains as many correlation peaks. Each path of the channel can be modelled 
by a complex multiplicative coefficient and a delay. The signals being propagated 
along the different paths can be aligned and combined by means of complex coefficients 
which are conjugates of the path coefficients, thus effecting a filtering adapted to the 
transmission channel. In order to simplify the terminology, the general expression 
"filtering adapted to the user tT will encompass both the filtering operation adapted to 
the signature of the user k and the filtering operation adapted to the transm.ssion 
channel. 

To combat the interference between signals destined for (the downlink) or coming 
from (the upUnk) the different users, multiuser detection methods have been proposed, 
and notably iterative detection methods such as those known as PIC (Parallel 
Interference Cancellation) and SIC {Serial Interference Cancellation). They are based 
on the iteration of an interference elimination cycle including the estimation of the 
symbols transmitted, the evaluation of the interferences and their subtraction from the 
signals received. Although of high performance, these methods are not optimal s.nce 
they do not provide an estimation in the sense of the maximum likelihood of the 
symbols transmitted by the different users. 

A method of multiuser detection with maximum likelihood inspired by the Viteib. 
algorithm was proposed by S. Verdu in an article entitled "Minimum probability of 
error for asynchronous Gaussian multiple access channels", published m IEEE 
Transactions on Information Theory, pages 85-96, January 1986, but its complexrty is 
prohibitive since it varies exponentially with the number of users. 

More recently a method of multiuser detection with maximum likelihood using a 
representation by a lattice of points was proposed by L. Brunei et al, in an article 
entitled "Euclidian space lattice decoding for joint detection in CDMA system- 
published in Proceedings of ITW, page 129, June 1999. According to this method, a 
vector characteristic of the received signal representing a statistic sufficient for the 
maximum likelihood detection of the symbols transmitted by the different users is 
determined. It is shown under certain conditions that the characteristic vector can be 
represented as the point in a lattice disturbed by a noise. The detection then consists of 
seeking the point in the lattice closest to the point corresponding to the vector received. 
However, the dimension of the lattice to be used generally being 2*K where K ,s the 
number of users, the number of points to be lested is still very high. To s.mphfy 
detection, it has been proposed to limit the search for the closest neighbour to the points 
in the lattice belonging to a sphere centred around the point received. Th.s s.mptified 
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detection method, referred to as the "method of detection by spheres", will be disclosed 
below: 

The context is a multiple access mobile telecommunication system with direct 
sequence spectrum spreading (DS-CDMA) comprising K users communicating 

synchronously with a base station. 

Let rf*(0 be the complex symbol sent by the user k at instant /. This symbol 

belongs to the modulation constellation A, used by the user k, which wUl also be 
referred to as the alphabet of symbols of the user k. Each user k transmits a block of 
N symbols with an amplitude of the signal o». The symbols are spread by a complex 
signature s t {() -s?(t) + jx s[ (0 with a duration equal to the symbol period T : 

s t (t)-0if tg[0,T] 

The * complex symbols d t {Q-d*(i)+jxdUi) transmitted at instant i are placed in 
a row vector of real values dj(r") defined as: 

The corresponding modulated signal is then, as a function of the time I : 

&jjifyz t d>(!)*(t-iT) < 2 > 

It is assumed that the channel is an ideal channel with white additive Gaussian noise. 
Let n-S,+tf be the signal received at time t and rj, a complex Gaussian noise of zero 
mean whose components have a variance N 0 ■ 

Let the row vector be y^W^^WI ^uch that jT«+-/*J<W is 

the complex output at instant / of the filter adapted to the user k : 

o 

The autocorrelation matrix of the spreading sequences will be denoted R(i). 

If the complex elements of (3) are decomposed into their real and imaginary parts, 
there is obtained; 
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Let A2~DUg(^^..,aM*) and Ra be the matrix of size 2K x IK such that: 



Equation (4) can then be put in matrix form: 

y2 (f)=da(0M2+n2(f) < 6 ) 

where M* is a real matrix of size 2K x2K defined by M2-A2R2 and where the noise 
vector n 2 (/>(^(/)^^%.^(iH(0) ^ N Q K 2 as its covariancc matrix. 

It will be demonstrated below that y 2 (/), as given by equation (6), can be 
represented as a point in a lattice A, of dimension 2*K, with a generator matrix Ma 
corrupted by a noise m . 

The term real lattice of points A of dimension wwill be used for any set of vectors 
of R" satisfying: 

x = b,v, + b,v 2 + ... + b.v, where b; 6 Z,Vi - I, ... ,k 
and where {vi,V2,,^vjc } is a base on R*. 

An example of a lattice of points of dimension 2 has been shown in Fig. 1. 

The points in the lattice form an additive abelian sub-group of R", and it is also 
the smallest sub-group of R r containing the vectors fovz^v*} and a Z-modulus of 
R r . These base vectors form the rows of the generator matrix G for the lattice. It is 
therefore possible to write x - bG where b~(&,...A)3* r . (7) 

The region delimited by the base vectors is referred to as the fundamental 
parailelotope and its volume, denoted vol(A) and det(A), is referred to as the 
fundamental volume. This fundamental volume is none other than the modulus of the 
vectorial product of the k base vectors and is therefore equal to |det(G)| where det 
designates the determinant. Though there are several possible choices for the generator 
matrix for the same lattice, on the other hand there is only one value for the fundamental 




volume. 
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The Voronoi region V or Dirichlet cell of a point x belonging to the lattice is all 
the points of R" closer to x than any other point in the lattice. The volume of this 
region is equal to the fundamental volume. 

The stacking radius p of the lattice is the radius of the largest sphere inscribed in 
the Voronoi region and the radius of coverage that of the smallest sphere circumscribed 
in this same region. The stacking radius is therefore the radius of the spheres whose 
stacking constitutes the lattice of points and the radius of coverage is that of the smallest 
spheres which, centred on the points in the lattice, make it possible to cover the entire 
space R' The density of the lattice is the ratio between the volume of the sphere of 
radius p and the fundamental volume. Finally, the coefficient of error (the kissing 
number) t(a) of the lattice is the number of spheres tangent to one and the same sphere 
in the stack or, in other words, the number of neighbours of a point in the lattice, 
situated at the minimum distance - 2p . 

Consider once again equation (6). The components of the vector d,(i) belong to a 
finite alphabet A of cardinal: 

The term A will be used for the constellation of the system (or simply constellation) 
as opposed to A 4 said to be modulation constellations. 

Assume for example that the components JP<i) and <//(») are PAM modulation 

symbols of order M: 

d?{i)E{-M+l,-M+3,..M-lM-'l} and (9) 
4(>)={-M + l -M^.M-lM-l) (10) 
If the transformation is effected: 

A W-J^tf >Jtf-0 «« d ^0>^0>^-O or a 8* n vectorially: 

d^HH>™> (n) 



where yju^M -1^-1,..^-!) 
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the components and d$) are elements of Z and consequently d'»(r) is a vector of 

In general terms, if there exists an affine transformation transforming the 
components dfff) and dl(i) into elements of Z, the vector t*(f) can be represented by a 

vector of Z 2K . 

In a similar manner, the corresponding transformation is effected on yi{i), that is 
to say: 

yi$H<y4>v«Mb) (12) 

By means of this transformation, which is assumed to be implicit hereinafter, the 
vector ch(0M2 then belongs to a lattice of points A 2 of dimension 2*K as defined by 
equation (7) with G-IVh. The vector y 2 (i) can then be considered to be a point in the 
lattice A 2 corrupted by a noise n 2 (/). 

If it is assumed that the components of the noise vector n 3 (/) are centred 
independent random Gaussian variables, the problem of detection in the sense of the 
maximum likelihood of the symbols transmitted by the different users amounts to a 
search for the point Zj in the lattice A 3 such that its distance to y,(/) is at a minimum. 

In reality, the components of the noise vector n,(i) are correlated and the 
covariance matrix of n 2 (i) is tf 0 R 2 . 

In order to boil down to the decorrelated case it is necessary, prior to the decoding, 
to effect an operation of whitening of the noise. 

The matrix R being her mi ti an, the autocorrelation matrix R 2 is symmetrical 
defined positive and can therefore be the subject of a Cholesky factorisation: 

R 2 -W 2 W 2 T ( l3 ) 

where W2 is an inferior triangular matrix of size 2£*2X 

A whitened observation vector is defined: yzt/^yz^^ 7 ^ 1 (14) 
and a new lattice of points Q 2 consisting of vectors of components 
(^}%(t%^t(ilxi(fb with ^(rJ-x^Owr 1 where 12(1) is a vector of components 
(xf(iixf^.,x{(f)jct(r)) belonging to A*. The lattice Q, has an inferior triangular real 
matrix as its generator matrix A 2 W 2 . 

It can easily be shown that, after whitening, the covariance matrix of the filtered 
noise n 2 (/)W 2 r " 1 is equal to MJ=r where l 2K is the identity matrix of dimension IK. The 
detection therefore comprises a first step of whitening the observation vector followed 
by a step of seeking the closest neighbour within the lattice of points Q 2 . 
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In order to reduce the number of points to be tested, as illustrated in Fig. 1, it is 
possible to limit the search to a sphere centred around the point yi . In practice, the 
choice of the radius of the sphere results from a compromise: it must not be too large in 
order not to lead to an excessively high number of points and sufficiently large to 
include at least the closest neighbour. 

Fig. 2 depicts schematically a multiuser detection device using a method of 
detection by spheres. The received signal n is filtered by a battery of filters adapted to 

each of the users, 210 ,210k. The real and imaginary components of the observation 

vector y,(/) output from the adapted filters are transmitted to a matrix calculation unit 
performing the spectral whitening operation according to equation (14). The real and 
imaginary components of the whitened vector yi(i) are then transmitted to a unit for 
detection by spheres seeking the closest neighbour of the point received within the 
lattice Q 2 of dimension 2*K. The coordinates of the closest neighbour directly give the 
real and imaginary components of the estimated symbols dk{i) for the different users. 

The step of searching for the closest neighbour is greedy in calculation time, 
which may prove very disadvantageous when the number of users is high. 

The aim of the present invention is to propose a simplification of the method of 
detection by spheres. 

To this end, the invention is defined by a method of detecting a plurality of 
symbols (d k (i)) transmitted by or for a plurality K of users, each symbol belonging to a 
modulation constellation and being the subject of a spectral spreading by means of a 
spreading sequence, the said method comprising a filtering step adapted for supplying 
a real vector (z) characteristic of the said received signal, at least the closest neighbour 
to the said vector being sought within a lattice of points (E) generated by the said 
modulation constellations, the said search being limited to candidate vectors (x), each 
of whose components (A) has a value lying within a search interval defined by a lower 
bound and an upper bound (B? ), the said bounds being chosen so that each of the 
said intervals comprises only values of components of points situated within a sphere 
of predetermined radius (JC ) and belonging to a modulation constellation. 

Advantageously, for a component with a given index i, each pair of lower (2F ) 
and upper {Bt ) bounds is calculated from a quantity T t characteristic of the width of 
the search interval for the said component. This characteristic quantity is determined 
by recurrence on the index /: the characteristic quantity (T t ) for a given index is 
determined from the quantity with the previous index (7; + ,) and the value of a 
component chosen in the search interval ([BT.i&n]) relating to the component 

with the previous index. 
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For a candidate vector component with a given index (2k\ the lower bound of 
the search interval will be advantageously chosen so as to be greater than the smallest 
integer (A/i) corresponding to a symbol of the modulation constellation (A k ) and the 
upper bound of the search interval will be chosen so as to be less than the greatest 
integer ( Mi* ) corresponding to a symbol of the modulation constellation. 

For a first vector candidate component with a given index (2*-l), the lower 
bound of the search interval will be chosen so as to be greater than a first integer 
( Mu-i ) and the upper bound of the search interval will be chosen so as to be less than 
a second integer (A/u-i ), the said first and second integers being determined from the 
value of a second component of the vector candidate such that the first and second 
components relate to one and the same user. 

Advantageously, the said first and second integers are respectively determined as 
the smallest value and the largest value of the first component such that the complex 
number defined by the value of the said first component and the value of the said 
second component is a symbol of the modulation constellation of the said user. 

The search for the closest neighbour is preferably carried out by scanning at the 
said search intervals one by one for the different component indices (/) and by 
choosing a component value (*,) in each of the said intervals, the bounds of each 
interval being determined according to the width of the interval with the previous 
index (/+ 1 ) and the component value (b H} ) chosen in this same interval. 

According to one embodiment, if, during the said search, the norm (|wfl) of a 

candidate vector is less than the radius of the said sphere, the radius is updated to the 
value of the said norm. 

The characteristics of the invention mentioned above, as well as others, will 
emerge more clearly from a reading of the following description given in relation to the 
accompanying figures. 

Consider once again a DS-CDMA telecommunication system with K users 
functioning in synchronous mode. As seen above, the detection of the symbols 
transmitted by the different users in the sense of the maximum likelihood can be 
reduced to a search amongst the points in a lattice (Q*)for the closest neighbour of the 
point corresponding to the received signal. 

In the case where the spreading sequences are real or more generally real multiples 
of the same complex number, it can be shown that the search can be carried out in a 
lattice with a dimension reduced to K. This is because, the imaginary terms of the 
matrix R 2 and consequently of the matrix Ma being zero, it is possible to boil down to a 
lattice of real points A of dimension K and of generator matrix M : 
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y*(/>d*(0M(/)+n*(/) 0 5 ) 

y'(j>d'(r)M(/)+i»J(f) ( 16 ) 

where y*<0,d*(0>tt*(0 (or respectively y^O^^/),!!^/)) are the vectors consisting of the 

reaJ parts (or respectively of the imaginary parts) of the components of yif^iOMO ; 

r 

M(f>AR(0 where R(r) is the matrix consisting of the coefficients /fa -p<J)&(f)df and 

A is the vector of the amplitudes of the K users. The observation vectors y*(i) and 
y ; (/) belong to R r . After any transformation according to an equation of the same type 
as (12), the vectors y*(i) and y'(f) can be considered to be points in a lattice A of 

generator matrix M(i) corrupted by noise. 

It can easily be shown that the noise vectors n*(/) and both have the 

covariance matrix Wo.R(f) . R being a symmetrical matrix defined positive, it can be 
factorised according to a Cholesky decomposition: R - W T where W is an inferior 
triangular real matrix of size KtlK. In order to decorrelate the noise components, the 
real observation vectors y R (i) and y'(0 are first of all subjected to a whitening 
operation: 

rowtiw*' 1 (17) 
?wr' < l8 > 

Secondly, the closest neighbours of the vectors y*(i) and y J (j) belonging to the 
lattice of points ft formed by the vectors x^)=x(j)W r " 1 , where x(0 belongs to A, are 
sought. It should be noted that the generator matrix of the lattice O is equal to AW, an 
inferior triangular real matrix. In addition, it can easily be shown that, after whitening, 
the covariance matrices of the filtered noises n*(j)W r ~ l and n'(f)w r ~ are both equal to 
Nolx where I* is the identity matrix of dimension K. 

When the symbols of or for the users are transmitted asynchronously, the 
modelling of the system is more complex since it is necessary to take account of the fact 
that a symbol of a user can interfere with two or even several consecutive symbols of 
another user. It can be shown in this case that it is possible to boil down to a search for 
the closest neighbour within a lattice of dimension 2*K % (K' in the case of real 
signatures) with K*>K where K 1 is a function of the number of symbols not yet 
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estimated which may interfere with each other. The detection is however not optimum 
in the sense of the maximum likelihood. 

In all cases, the problem is to determine the point x in a lattice S of dimension k 
closest to the received whitened vector y , which amounts to minimising the metric 

where y-x+Tj, ri*4n±-Sl*) is the noise vector and x-(»,...,Xr) is a point 

belonging to the lattice. The noise vector rj has independent real components in a 

Gaussian distribution of zero mean. 

It should be noted that the vector y(/) does not need to be whitened if use is made 

of a metric based on the covariancc matrix: 

/«(y/x)-(y-x)R- 1 (y-x) r (20) 

Hereinafter, for reasons of simplification, the observation vector, whhened (y(0) 
or not (y(0) will ^ termed z and the metric acting in equation (19) or (20) will be 
termed ||| . 

The points in the lattice S can be described by the vectors x=bG where 
h^{bi^.Jbk) has components b t belonging to the ring of the integers Z and where G is 
the generator matrix of the lattice. The rows of the matrix G are denoted {vi,w ■» v*}. 

By definition these vectors form a base of the lattice. 

The set of transmitted symbols is limited to an alphabet of finite size AkCZt 

referred to as a constellation. This constellation is determined by the modulation 
constellations used by (or for) the k users and the cardinal of the alphabet Ar is the 
product of the cardinals of the different modulation alphabets. It is assumed that the 
complex points of each of these constellations have real values and evenly distributed 
imaginary values. 

As has been seen, an exhaustive decoding would require a search for the closest 
neighbour throughout A*. The decoder advantageously restricts its calculation to the 
points which are situated within an area of the constellation situated around the received 
point, preferably inside a sphere of given radius -JC centred on the received point as 
depicted in Fig. 1. Only the points in the lattice situated at a quadratic distance less than 
C from the received point are therefore considered for the minimisation of the metric 
(19). 

In practice, the decoder effects the following minimisation: 
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To do this, the decoder seeks the smallest vector w in the translated set z-E . 
The vectors z and w can be expressed as: 



It is important to note that p and £ are real vectors. As w~z-x, where x 
belongs to the lattice E, this gives the equation f , - p, - for /-l with w« y§<v, 



The vector w is a point in the lattice whose coordinates £ are expressed in the 
translated reference frame centred on the received point. The vector w belongs to a 
sphere of quadratic radius C centred at 0 if: 



Ln the new system of coordinates defined by 5, the sphere of quadratic radius C 
centred at y is therefore transformed into an ellipsoid centred on the origin. The 
Cholesky factorisation of the Gram matrix T=GG r gives T-AA T , where A is an 
inferior triangular matrix of elements 6V . 

Il should be noted that, if the vector y has been whitened, it is not necessary to 
effect this factorisation since the generator matrix for the latuce is already triangular and 
inferior. 

However, where prior whitening has not been carried out, and therefore where 
Cholesky decomposition is necessary: 



z=gG 
w-gG 



with e=o>->/*) 
with (-&.--£) 
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By putting 



5 U 2 for i 
^ for j 



1 , ... , X, 



1, 



1 rz j + 1, ... ,X 



(25) 



there is obtained 



200 2-2 1 7 8 6 9 . 



(26) 



Being concerned first of all with the range of possible variations oF , and then 
by adding the components one by one, the following k inequalities are obtained, which 
define all the points within the ellipsoid: 



(27) 



It can be shown that the inequalities (27) make it necessary for the integer components 
of b to satisfy: 



(28) 



where [x] is the smallest integer greater than the real number x and [x] is the largest 

integer smaller than the real number x . 

The decoder has k internal counters, namely one counter per dimension, each 
counter counting between a lower and upper bound as indicated in (28), given that each 
counter is associated with a particular pair of bounds. In practice these bounds can be 
updated recursively. We put: 

SSfoi* ■ + Jj^ (29) 



(30) 
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T^-Tt-qfa-bJ (3D 

whh r*-c 

Using equations (29) to (31), the range of variation of each component b t is 
determined recursively, commencing with the component k : 

U*h±ti (32) 



with Z7« 



and Zi - 



\qtt 



(33) 



The ranges of variation defined by (33) are advantageously restricted so as not to 
unnecessarily test points which are situated outside the constellation. It should be stated 
that each user k uses a constellation of symbols A k , each symbol of the constellation 
being a complex number whose real and imaginary parts (possibly after affine 
transformation) are elements of Z. the general case is first of all taken where the lattice 
is of dimension 2K t where AT is the number of users. For each user k 9 the dimensions 2k 
and 2k-\ bearing the complex symbol transmitted by or for the user k are considered. 
As indicated in Fig. 3, the constellation, or equivalently the modulation constellation of 
the user A, is first of all projected onto the dimension 2k This projection defines an 
interval [v/u^A/i]. The search interval [Bik^u] is then defined where: 

B' 3k - Upp(L- 3k Mh ) *U « Lo»(U 3k MJ k ) (34) 
and an integer b 2k is chosen in this interval. The component b^ then defines an interval 
[MZk-iM**-*H for the component b 2kA as indicated in Fig. 3. A search interval 
[J92*-vEz*-i] is then defined by: 

B~ 2k _ 3 - Vpp(L- 2k _ } M; k .j ) and B] k _ s - Low(U 7k ^ ) (34') 

By proceeding in this way it is ensured that the closest neighbour is sought only 
amongst candidates which are both situated in the search sphere and points in the 
constellation. 

The situation is slightly different when real spreading sequences are used and 
therefore when the lattice is of dimension K (where K is the number of users). In this 
case, the real and imaginary values of the symbols are the subject of parallel searches in 
a lattice (A, Q) of dimension K. The calculation of the search bounds (and the choice of 
the values of the components) relates by turns to the real part and the imaginary part, in 
an interleaved fashion. 

In the simple example where the modulation employed on a dimension i is a PAM 
modulation of order M t the integer coordinate b t of the point sought must be between 0 
and M-l. The search bounds are adjusted by B: ^UppfL'fl) and 
B* =Lvw(L\M -1) so that the counter associated with the component b t does not 
travel over points which are situated outside the search sphere or outside the 
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constellation. This adjustment of the search bounds considerably accelerates the 
algorithm for detection by spheres. 

In addition, ihe search within the sphere can be accelerated still further by 
updating the radius VC with the last calculated Euclidian norm(|wj . 

The initial value of the search radius Vc must be chosen in an appropriate 
manner. This is because the number of points in the lattice situated within the decoding 
sphere increases with C. This is why the choice of a large value of C disadvantages 
the decoding algorithm whilst the search sphere may be empty if C is too low. 

So as to be sure that the decoder finds at least one point in the lattice, a search 
radius greater than the radius of coverage of the lattice is advantageously chosen. It can 
for example be taken so as to be equal to the upper Rogers bound: 

^-(iriog.c+rioglog^+S^^l^ (35) 

where Vk is the volume of a sphere of unity radius in the real space R* . 

Fig. 3 depicts a flow diagram for the search for the closest neighbour used in the 
method of detection by spheres according to the invention. 

The vector 2, the matrix G and the initial value C are the parameters transmitted to 
the search method. 

It will be assumed that the situation prevails where k~2K, the simplified case 
jr=X: being derived therefrom without difficulty. 

In a first step 301, the Cholesky factorisation of the Gram matrix T-GG 3 " is 
carried out if the vector z has not resulted from a whitening (otherwise the matrix G is, 
as has been seen, already triangular and inferior). The coefficients ^ are calculated 
according to equation (25) and the values pt are calculated by the equation: g-zG" . 

At step 302 the quadratic value of the smallest current distance is initialised to C: 
££»-C and the recursive calculation of the search bounds is initialised by: 7i«=C , 
£-/M-l v ../r • The values of bounds Mi ,Mt are also determined with i-2k 9 k=\,„,K 
and next K tables are constructed giving the values Kdi ,Mt with i=2k-J y k=]...yK 
corresponding to the different possible values b^ 

At step 303 the index of the search dimension is initialised, that is to say f-r . 

At step 304 the values U ,B are calculated by means of equation (33). ET ,Bt 
are determined by means of equation (34) if / is even. If i is odd, M" \Mt are 
determined from the table (/+l)/2 and the value of b m . The hounds are next 

determined according to equation (34'). In addition, b f is initialised to the value £"-1 . 

At step 305 the value of b, is incremented by 1 : 6-6+1 . 
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Step 306 consists of a test: £»2T ? 

If the test at 306 is negative, step 307 is passed to and h is tested whether i>\. 
If the test at 307 is positive, the search bounds are calculated step by step. At step 
310, r^-TI-^te-Ai) 2 and £, - A - ^ calculated and then, at step 311. 

&-i-p-i+^9*-i£j is calculated. Next i is decremented at 312 before returning to step 

304 in order to calculate the search bounds for /-I . 

When the test at 307 is negative, step 308 is passed to in order to calculate the 
quadratic norm |]w|f of the vector w relating to the current point by 
|(wj| 2 =^-fi+^£i-&) 2 . Then it is tested whether ||w|f<i&.. If such is not the case 

305 is continued with. On the other hand, if the minimum distance is improved, step 
313 is passed to, at which the components b t of the vector x corresponding to this 
minimum are stored. The minimum quadratic distance is updated by DL=§wf and the 
size of the search sphere by K-||wf . The search procedure is recommenced by 
returning to 303, 

If the test at 306 is positive, that is to say if the upper search bound is reached for 
the current dimension, it is tested at 321 whether /-jr. If such is the case, the search is 
terminated at 320 and the last vector x stored is the closest neighbour. On the other 
hand, if the test is negative, the search continues by incrementing / at step 322 and 
returning to step 305. 

4. Brief Description of Drawings 

Fig. 1 depicts a lattice of points useful to the detection method employed in the 
receiver illustrated in Fig. 2; 

Fig. 2 depicts schematically the structure of a multiuser DS-CDMA receiver using 
a method of detection by spheres; 

Fig. 3 depicts a flow diagram for the search for the closest neighbour used in the 
method of detection by spheres according to the invention; 

Fig. 4 depicts an example of a user modulation constellation. 
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FIG.3 
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FIG. 4 



1. Abstract 

Method of detecting a plurality of symbols (d k (0) transmitted by or for a 
plurality K of users, each symbol belonging to a modulation constellation and being 
the subject of a spectral spreading by means of a spreading sequence, the said method 
comprising a Filtering step adapted for supplying a real vector (z) characteristic of the 
said received signal, at least the closest neighbour of the said vector being sought 
within a lattice of points (S) generated by the said modulation constellations, the said 
search being limited to candidate vectors (x), each of whose components (£,) has a 
value lying within a search interval defined by a lower bound (2JT ) and an upper bound 
( R+ ), the said bounds being chosen so that each of the said intervals comprises only 
component values of points situated within a sphere with a predetermined radius (y/C ) 
and belonging to a modulation constellation. 

2. Representative Drawings 

Fig. 3 



